Spin-flop transition in antiferromagnetic multilayers 
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A comprehensive theoretical investigation on the field-driven reorientation transitions in uniaxial 
multilayers with antiferromagnetic coupling is presented. It is based on a complete survey of the 
one-dimensional solutions for the basic phenomenological (micromagnetic) model that describes the 
magnetic properties of finite stacks made from ferromagnetic layers coupled antiferromagnetically 
through spacer layers. The general structure of the phase diagrams is analysed. At a high ra- 
tio of uniaxial anisotropy to antiferromagnetic interlayer exchange, only a succession of coUinear 
magnetic states is possible. With increasing field first-order (metamagnetic) transitions occur from 
the antiferromagnetic ground-state to a set of degenerate ferrimagnetic states and to the saturated 
ferromagnetic state. At low anisotropies, a first-order transition from the antiferromagnetic ground- 
state to an inhomogeneous spin-flop state occurs. Between these two regions, transitional magnetic 
phases occupy the range of intermediate anisotropies. Detailed and quantitative phase diagrams 
are given for the basic model of antiferromagnetic multilayer systems with A = 2 to 16 layers. The 
connection of the phase diagrams with the spin-reorientation transitions in bulk antiferromagnets 
is discussed. The limits of low anisotropy and large numbers of layers are analysed by two different 
representations of the magnetic energy, namely, in terms of finite chains of staggered vectors and in 
a general continuum form. It is shown that the phenomena widely described as "surface spin-flop" 
are driven only by the cut exchange interactions and the non-compensated magnetic moment at the 
surface layers of a stacked antiferromagnetic system. 

PACS numbers: 75.70.-i, 75.50.Ee, 75.10.-b 75.30.Kz 



I. INTRODUCTION 

Since the discovery of antiferromagnetic interlayer 
exchangei and the giant-magnetoresistancaS in magnetic 
superlattices, such structures have become important 
components in magneto-electronic devices. Research on 
these coupled multilayer systems is mainly driven by 
applications in magnetic storage technologies and the 
emerging spintronics.^ Specific structures with antiferro- 
magnetic coupling are now considered as promising stor- 
age medial^ It is clear, that applications necessitate a 
thorough control and understanding of their magnetic 
properties. On the other hand, such synthethic anti- 
ferromagnetic structures are ideal experimental models 
for studies of magnetic states and magnetization pro- 
cesses of antiferromagnets in confining geometries 
Two ferromagnetic layers coupled antiferromagnetically 
through a spacer, as the simplest of these systems, show 
properties which are formally described by the same phe- 
nomenological theory as a two-sublattice bulk antiferro- 
magnet. However, the magnetic states, domain struc- 
tures, and magnetization processes even of such two- 
layer systems display a bewildering variability and are far 
from understood in detail.^ Modern experimental meth- 
ods now allow imaging of magnetic states and domains 
in multilayer systems with resolution into the depths of 
multilayer stacksiiSiii Therefore, detailed studies of such 
structures have become feasible. 

Theoretical models to describe the magnetic states 
of finite antiferromagnetic superlattices have revealed 
various surface effects, rich phase diagrams, and com- 
plex magnetization processes. For antiferromagnetic lay- 



ers, there are many other effects. Surface-induced in- 
teractions, exchange couplings of antiferromagnets to 
other magnetic systems, in particular exchange bias in 
antiferromagnetic-ferromagnetic bilayer systems^^ add to 
the multitude of possible magnetic states in antiferro- 
magnetic layers^ However, the difficulty to probe and 
image magnetic structures in antiferromagnetic materials 
impedes the progress of our understanding on the anti- 
ferromagnetic side in such layered systems. Hence, the 
finite antiferromagnetic superlattices are a suitably sim- 
ple system which may promote a better understanding of 
surface related effects in antiferromagnetic layered sys- 
tems generally. It is important to stress here, that sur- 
face effects in antiferromagnets have a different nature 
than in a ferromagnetic system. The cut exchange bonds 
at a (partially) uncompensated surface of an antiferro- 
magnet causes a particular disbalance of magnetic forces 
which can never be understood as a small surface-effect. 
In contrast, surface-effects in ferromagnets are related to 
spin-orbit effects which usually are weak in comparison 
to the exchange. 

A stack of magnetic layers with antiferromagnetic cou- 
plings provides the basic model for cut exchange bonds 
at a fully uncompensated surface. The study of the mag- 
netic states and transitions for such systems in exter- 
nal fields has a long history. In 1968 Mills proposed 
that, at the surface of a uniaxial antiferromagnet, a first- 
order transition should occur in fields below the common 
"bulk" spin-flop (SF). This transition from the antiferro- 
magnetic to a "surface spin-flop" state should result in 
flopping a few layers of spins near the surface, i.e. they 
would turn by nearly 90 degreeii^ Further theoretical in- 
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FIG. 1: (Color online) Sketch of an antiferromagnetically 
coupled multilayer corresponding one-dimensional spin chain. 
The "exchange springs" are cut at the ends of the finite chain. 



vestigations have improved the mathematical analysis of 
this reorientation effecti^iiSiiLiSiiS however, direct exper- 
imental observations at surfaces of crystalline antiferro- 
magnetic materials failed, e.g., for the classical uniaxial 
antiferromagnet MnF2, (see bibliography and discussion 
in 113). 

As was mentioned above the basic model for two an- 
tiferromagnetically coupled layers is equivalent to the 
classical mean-field description of two-sublattice bulk 
antiferromagnetSsSiiSS, These systems compose two large 
groups: antiferromagnets with weak anisotropyS^ and 
strongly anisotropic uniaxial crystals that are commonly 
called metamagnets.^^ The metamagnetic phase transi- 
tion between the antiferromagnetic and ferromagnetic 
phases has been observed and investigated in many an- 
tiferromagnetic bulk systemsi^ Easy-axis antiferromag- 
nets with weak anisotropy also compose a large group 
of magnetically ordered crystals. For bulk antiferro- 
magnets, the spin-flop transition has been predicted by 
Neel2^ and later was observed experimentally in CuCl2- 
2H 20i24 In the next fifty years, spin flop transitions have 
been discovered and carefully studied in ma ny classe s of 
antiferromagnets (see bibhography in Refs. |2a25l26j ). 

The interest in spin-flops revived with the synthesis of 
magnetic multilayer stacks with indirect antiferromag- 
netic exchange coupling through spacers,^ These artificial 
antiferromagnetic layers, with few magnetic units as 
macroscopic spins instead of atomic spins (Fig. ^ 
offer the possibility to study field-driven reorientation 
transitions with unusually low exchange compared to 
anisotropics.^ Experimental investigations in Fe/Cr(211) 
antiferromagnetic superlattices^ seemed to confirm the 
scenario of the surface spin-flop transition introduced 
m Refs. Pll^ . The observation was supported by 
numerical investigation of an elementary micromagnetic 
modeliSiS& But, investigations on the magnetism of 
such antiferromagnetic superlattices and thin films 
did not produce a consistent understanding of the 



occurrence and nature of "spin-flop" transitions or other 
reorientation transitions in layered antiferromagnetic 
structures. Recent experiments demonstrate complex 
behaviour and different scenarios for the evolution of the 

magnetic StateS,l°i"il^1^7l^^1^9,30,31,32,33,34,35,36,37,38,39 

The problem of the magnetic states in antiferro- 
magnetic superlattices with uniaxial anisotropy 
is strongly related to the long standing prob- 
lem of a "surface spin-flop" discussed in various 
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oretical investigations are mostly based on numerical 
calculations within a certain model, here called Mills 
modeZi'^^i^^i^^i^^i^^i'^'^ Further theoretical works related 
these finite or semi-infinite antiferromagnetic chain mod- 
els to systems like the Frenkel-Kontorova model42iS 
These studies have led to controversial results and gener- 
ated a long-drawn discussion about the physical nature 
of the reorientation transformations in the antiferro- 
magnetic superlatticesi'^^i'^'-'i^^i^^i'^^i^^ Only few attempts 
have been made to gain a complete understanding of 
the ground-state structure of the basic one-dimensional 
models for these superlattices42i^2iS Thus, the basic 
questions, whether, when, and how a surface spin-flop 
occurs, were unresolved. 

A full set of solutions for the model under discussion 
was recently obtained by usiS These outwardly simple 
systems with few degrees of freedom own rich phase dia- 
grams because of the competition between internal stiff- 
ness and anisotropy in conjunctions with restricted di- 
mensionality. Present paper presents an exten ded ac - 
count and an analysis of the solutions from Refs. |53l54| . 
We explain the physical mechanism responsible for the 
formation of the main magnetic states in antiferromag- 
netic superlattices. We derive a clear and simple picture 
of the phenomena which have been discussed as "surface 
spin-flop" . We demonstrate the connections with bulk 
antiferromagnetic systems and other classes of magnetic 
nanostructures .^^1^^ 

The structure of the paper is as follows: The phe- 
nomenological model and its variants are introduced in 
Sec. m The analysis of possible magnetic states in the 
system starts from two limiting cases of low and high uni- 
axial anisotropy. The full solution for the generic highly 
symmetric Mills model in applied fields along the easy 
axis are presented in Sec. IIIII and the general structure 
of the phase diagrams are presented. In this section, 
mainly analytical and numerical examples are employed 
to explain these solutions. Generalizations are briefly in- 
dicated. The methods can be extended to solve any other 
more general model for an antiferromagnetic multilayer 
system. In Sec. IIVI we exploit the fact that the models 
for low-anisotropy antiferromagnetic superlattices can be 
reduced to equivalent models for a chain of exchange cou- 
pled two-sublattice antiferromagnets. This dimerization 
transformation reveals the physical mechanisms ruling 
the magnetic states and reorientation transitions. Nu- 
merical results for the evolution of inhomogeneous spin- 
flop phases, magnetization curves, and phase diagrams 
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are presented along with this discussion. In particular 
the limit of large numbers of layers and the emergence of 
the bulk spin-flop for such finite systems are discussed. 
Then, the continuum representation of the general mod- 
els is presented. This offers a different point of view for 
the weak anisotropy. This continuum approximation is 
applicable to any weakly anisotropic system. It allows 
to derive the structure of the inhomogeneous spin-flop 
phase for arbitrary models, in particular model where the 
magnetic moments at the surfaces are partially compen- 
sated. Hence, the succession of field-driven phase tran- 
sitions between the antiferromagnetic state towards the 
saturated ferromagnetic state via this spin-flop phase can 
be completely analysed. In Sec. IV Al the general picture 
of magnetic states in antiferromagnetic nanostructures 
and some recent experimental results are discussed in the 
context of the new results. 



II. MODEL 

A. The micromagnetic energy. Mills model 

Let us consider a stack of N ferromagnetic plates infi- 
nite in X- and y-directions and with finite thickness along 
the z-axis. The magnetization of each plate is mi, and 
they are antiferromagnetically coupled through spacers. 
Replacing this system by a chain of single-domain par- 
ticles with spontaneous magnetization TOq'^* — |mi|, we 
may describe the magnetic configuration by the variables 
Si = mi/mg*\ i.e., by the set of unity vectors along the 
magnetization of the ith layer. We assume that the fer- 
romagnetic layers have a uniaxial magnetic anisotropy 
with a common easy axis. The phenomenological energy 
of this system can be written as 

N-l N 

©Af = ^ ^Ji Si • Si+i -I- Ji (Si • Si+i)^ - H • ^ Ci Si 

1=1 1=1 
^ N N~l 



Here, Q — m^o' /tjiq designate deviations of the mag- 
netization in the i-th layer from the average value mo. 
Ji and Ji are constants of bilinear and biquadratic ex- 
change interactions, respectively. The unity vector n 
points along the uniaxial anisotropy direction; Ki and 
JsTj' are constants of the in-plane and inter-plane uniaxial 
anisotropy. Finally, ea{va.i) collects higher-order uniaxial 
and in-plane magnetic anisotropy contributions, e.g., in- 
trinsic cubic (magnetocrystalline) anisotropy in systems 
like Fe and Ni layers. 

The functional ^ generalizes similar models con- 
sidered earlier in a number of studies on magnetic 
states in antiferromagnetic multilayers with uniaxial 
anisotropyi^i^°i'^°i^^i'*^i^'^i'^^i^^i'^''i'^^i^^ In a recent paper on 
antiferromagnetic superlattice with higher (tetragonal) 



symmetry, Ref. [Sa ] , we have discussed the general model 
and the justification of the approach used for this class 
of systems. Therefore, the reader is referred to^S- for a 
detailed discussion and further references. Eq. com- 
prises the magnetic energies due to the main physical 
mechanisms, which are present in magnetic multilayers 
with indirect exchange through spacer layers (see Fig.^l. 
The ferromagnetic layers can be considered as homoge- 
neously magnetized blocks with constant values of the 
magnetic interactions. This assumption relies on the fact 
that in ferromagnetic nanolayers the intrinsic (direct) ex- 
change coupling are usually very strong. Thus, they play 
the dominating role for the magnetic order within the 
layers which react also very stiffly on all external and 
induced magnetic forces. Thus, the internal magnetic 
structure of an individual layer experiences little change 
under infiuence of the induced magnetic forces at the sur- 
faces and interfaces, and the reorientation of other layers 
in the stack. This hypothesis has been justified by suc- 
cessful applications of such models to describe magneti- 
zation processes in layered ferro- and antiferromagnetic 
nanostructures in different classesi^i-^-^i-^'^i^^i^'^'i^^i^^i^" 

Antiferromagnetic multilayers with uniaxial anisotropy 
and in applied field along the axis n show the strongest 
reorientation effects. Here, we address the overall mag- 
netic properties of such uniaxial systems. Their be- 
haviour may be analysed by considering only terms with 
bilinear exchange J, uniaxial intralayer anisotropy K 
and an external field H. The form of the energy in 
Eq. JQ) considers additional terms, which are known 
to play a role in antiferromagnetically coupled multi- 
layers. In particular, strong biquadratic exchange has 
been revealed in a number of antiferromagnetically cou- 
pled multilayers i^iiSiSiiSS In Ref. [s^, we have stud- 
ied the related phenomenological model for multilayers 
with zero and with four-fold (tetragonal) anisotropics and 
we have discussed the relevance of biquadratic exchange 
(J 7^ 0), which plays essentially only a quantitative role 
for the competition between the various possible mag- 
netic states as long the antiferromagnetic ground-state 
remains coUinear. 

The uniaxial anisotropics may be intrinsic to the mag- 
netic material of the film or induced by surface effects. 
Thus, the uniaxial anisotropics can be strongly changed 
with respect to bulk systems, and their strengths can be 
controlled in film system within wide margins. Combina- 
tion with intrinsic and induced fourth-order anisotropy, 
as considered in Cq in Eq. (^), may lead to peculiari- 
ties of magnetic properties, see Refs. p8ll55l |. but here 
we will disregard these contributions. Finally, strong de- 
magnetization effects in antiferromagnetic superlattices 
with perpendicular anisotropy are responsible for com- 
plex evolution of multidomain states and specific magne- 
tization processeSiiiSi^SiSiSi To investigate these effects 
the stray-field energy must be included in Eq. , and the 
corresponding magnetostatic problem has to be solved. 

The magnetic superlattices with antiferromagnetic 
coupling can be separated into two classes: non- 
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compensated systems with a net magnetization and those 
with fully compensated magnetization. In the former 
case, the net magnetic moment strongly determines their 
global magnetic properties. In many physical aspect their 
properties are similar to those of bulk ferrimagnets. The 
main subject of this paper are the multilayers with fully 
compensated magnetization, i.e., multilayers with even N 
and equal magnetization in all layers are similar to bulk 
collinear antiferromagnets. For simplicity, we assume 

= 1 in all layers i = 1 . . . N. At the end of the paper 
we also consider effects imposed by a partial compen- 
sation of the magnetization in the cndmost layers, i.e., 
deviations of Ci and from unity. 

For the reorientation effects in the antiferromagneti- 
cally coupled chain (Fig. the effect due to the cut 
exchange bonds at the surface dominates. The last mo- 
ments in the chain are coupled only to one neighbour 
while "internal" moments interact with two. Thus, the 
moments at the surface experience a weakened exchange 
stiffness and are more susceptible to the reorienting in- 
fluence of an applied field. Due to this cut exchange at 
the surfaces the boundary moments oriented against the 
external field turn into the field direction in lower fields 
than internal moments. The simplified version of the 
model (P) with equal constants Ji — J, Ki — K, and 
Ji = K[ = Ca — Q for an applied field in direction of the 
easy axes, H || n, describes the effect of the cut bonds as 
the sole surface-imposed factor. It allows to investigate 
this surface effect separately from other interactions. 

Usually the magnetization of the layers is confined to a 
certain plane. For many multilayers, this is the film plane 
owing to demagnetization. For this case the deviations of 
mi from the anisotropy axis can be described by angles 
Oi, and the energy Q in the equal constant model is 
reduced to the following form 



Af-l 

$Ar = J > cos( 

i=l 

N 



z=l 

A 



(2) 



K 



N 



cos^ Oi 



i=l 



For J > and K > 0, the Eq. Q de- 
scribes antiferromagnetically coupled ferromagnetic lay- 
ers with easy-axis anisotropy in an external field 
along the anisotropy axis. Energy (O has been in- 
troduced by Mills for a semi-infinite chain [N = 
oo).^'* Later this model (called here Mills model) 
has been intensively studied for finite and infinite 

_A r5. 6. 15. 16. 17. 18. 19. 20.40. 41. 42. 43. 44. 45.46. 47.48.49. 50. 51 

been used as basic ansatz to analyze experimental results 
in antiferromagnetically coupled multilayersi^*ifi*^ 

Both the exchange interactions and uniaxial anisotropy 
in ||2J) have surface/interface induced nature, their values 
are very sensitive to many physical factors such as the 
structure of the spacers and substratesiSSiSi Hence, the 
materials parameters may vary from layer to layer in the 
stack. A generalization of the Mills model may include 



differing parameters for each layer 

Af-l 



(3) 



i=l 



N 



cos 



1 ^ 



cos 9i 



i=i 



This is a general model for antiferromagnetic multilay- 
ers. Even for a stack of identical nanolayers, the top and 
bottom layers still have a different "neighbourhood" than 
internal layers. To describe these effects we may intro- 
duce a modified Mills model with equal materials param- 
eters for all internal layers Ji = J {i = 2,3, . . . , N — 2), 
Ki — K {i — 2,3, . . . ,N — 1) and different values for the 
first and last layer (Ji Jn-i — Js, Ki — Kn — Ks, 

Cl = = Ca < 1) 



N 



Js [cos(6'i - 62) + cos{9n-i - 9n] 



(4) 



N-2 



-J ^ cos(6'j - 6*^+1) 



-H(s (cos 9i + cos 9n) — H ^ cos 9i 



i=2 



Ks (cos^ 6*1 + cos^ 9n) 



K 



N-l 



cos^ 9i 



i=2 



In all these cases, calculations of the magnetic states for 
the antiferromagnetic superlattices can be reduced to the 
minimization of the energy functions <I>Ar(0i, 92, ... , 9^)- 
In this paper, we study in detail solutions for the chains 
with equal constants as described by Mills model 
and we discuss the generalizations according to Eqs. 
and I©. 



B. General features of the solutions. Relation to 
bulk antiferromagnetism 

The antiferromagnetic multilayers with N = 2 
are of particular importance for investigations on 
surface/interface-induced interactionsiS2i2i422i2& In exper- 
imental works they are often referred as "trilayers" , we 
use here the term "two-layer systems". The energy Q 
for = 2 is the same function as the mean-field magnetic 
energy of a bulk two-sublattice antiferromagnet 



$2 = Jsi-S2~ H [(si • n) + (S2 • n)] 



K 

y 



(si • n)" + (S2 • n) 



(5) 



We briefly review the reorientation transitions in these 
two-layers systems to fix our notation and terminology. 
Following a general convention we introduce the linear 
combinations of the magnetization vectors S1.2 = M±L, 
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the total or net magnetization M and the staggered mag- 
netization L, that is also called vector of antiferromag- 
netic ordev^. The equations |si| — 1 lead to the con- 
straints -I- = 1 and (M • L) = 0. In a magnetic 
field along the easy axis the vectors H, Si and S2 lie in 
a fixed plane. The magnitude of the net magnetization 
M = |M| can be used as an internal parameter. After an 
independent minimization with respect to M, the energy 
and the net magnetization depend only on the orienta- 
tion of the staggered vector L and can be expressed as 
functions of the angle (j) between anisotropy axis n and 



M = 



{2J + K cos2<j)) 
2Hsin(j) 
(2 J + X cos 20) ' 



Kcos^cIj+J, (6) 
(7) 



Then, the phase diagram of the antiferromagnetic two- 
layer system is given by the potential energy for a system 
with one variable </> (Q with the control of materials pa- 
rameters K/J,H/J. This phase diagram is plotted in 
Fig. 121 The structure of the phase diagram is determined 
by the following characteristic fields 



i/tr = VK{2J-K), Hfm = J,Hf = (2 j - K) 



Haf = y/K{2J + K), HsF = (2 J - K) 



K 



2J + K 



(8) 



For K < J a first-order transition between the antifer- 
romagnetic (AF) phase with </() = 0, tt and the spin-flop 
phase (SF) with = ±7r/2 occurs at the spin- flop field 
HtY- The fields Haf, Hsf are stability limits of these two 
competing phases. The difference Haf ~ Hsf gives the 
width of the metastability region. Because the antiferro- 
magnetic state has zero magnetization the magnetization 
jump at the first-order transition, AM — mo AfsF(i?tr) = 
mo K/ (2 J — K) exactly in the field i?tr, equals the 
magnetization of the spin-flop phase. For increasing 
anisotropy K this magnetization jump Msf{Hsf) grad- 
ually increases from very small values M <C 1 for A' ^ J 
to the saturation value M — \ which marks the end point 
of the first-order transition line between the AF and SF 
phase that is reached at K = J . A continuous transi- 
tion from the SF phase to the ferromagnetic (F) phase 
occurs at iJp- This transition leading to an enforced field- 
polarized state is usually referred to as spin-flip transi- 
tion. 

For K > J, the SF phase does not arise as a sta- 
ble state; instead there is a direct first-order transition 
between the AF and F phase at Hi — J. Such tran- 
sitions in antiferromagnets are known as metamagnetic 
transition. For this high anisotropy region, K > J, the 
critical field Hf plays the role of the stability limit for the 
ferromagnetic phase (Fig. The metamagnetic transi- 
tion is characterized by a large jump of the magnetiza- 
tion AM (Hi) — mo and extremely broad metastability 
regions. 
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FIG. 2: (Color online) The phase diagram of the solutions 
for a two-layer system includes antiferromagnetic (AF), spin- 
flop(SF) and ferromagnetic (F) phases. The AF phase is sep- 
arated from the SF and F phases by the first-order transition 
lines (thick). The second-order "spin-fiip" line .ffp (dashed) 
separates SF and and F phases. Thin lines indicate stability 
limits of the corresponding phases. Inset: in the spin-flop field 
H — Hti the potential wells switch from (0, tt) (AF phase) to 
(-7r/2, Tv/2) (SF phase). 



For low-anisotropy systems in the region K <ti J, the 
energy (Q can be simplified 

$2(0)= (^ ^'^'^^" ^ cos 20, H„ = V2Jk. (9) 

In this limit the metastable region is restricted to a close 
vicinity of the spin-flop fleld: Haf ~ H^f ~ Hn ~ Hq. 
This simplified potential energy Eq. Q reveals the phys- 
ical mechanism of the spin-flop transition. At zero field 
the uniaxial anisotropy stabilizes the antiferromagnetic 
phase. The potential wells at = 0, tt corresponding to 
the two antiferromagnetic states are stable. An increas- 
ing applied field, H < Htr, gradually reduces the height 
of the potential barrier between the antiferromagnetic 
states. At the threshold spin-flop field Hu the stable po- 
tential wells switch into the configurations — —tt/2, 
= 7r/2 that correspond to the flopped states (Fig. 

Contrary to natural antiferromagnetic crystals which 
are described only by marginal parts of the K/J scale 
with either low anisotropy, K/J <C I, or very large 
anisotropy. A' 3> 1, antiferromagnets composed of mag- 
netic nanolayers can have arbitrary values of K/J. These 
artificial antiferromagnets cover the whole phase plane in 

Fig.ia 

In finite multilayers with N > 2 the cutting of the ex- 
change bonds at the surfaces (Fig. ^) causes a strong 
disbalance of the exchange interactions along the chain. 
This disbalance is the determining factor for the appear- 
ance of magnetic states in the system. The detailed 
analysis of the solutions for Mills model will be given 
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in the next section. Here, we summarize the physical 
mechanism ruhng the formation of magnetic configura- 
tions in simple terms. In the antiferromagnetic configu- 
ration the moments at one surface always point against 
an applied field. These moments can be reversed more 
easily than the internal moments because of the cut ex- 
change. Depending on the relative strength of the ex- 
change and anisotropy, this specific instability leads to 
different reoriented configurations. At very strong uniax- 
ial anisotropy (K 3> J) the exchange coupling between 
layers becomes negligible. In this case the reorientation 
of the magnetization in the endmost layer does not influ- 
ence magnetic states in the other layers. In an increasing 
field, a collinear spin configuration with an inverted end- 
most moment is reached. The corresponding ferrimag- 
netic (FM) phase becomes energetically stable at a field 
^FM = J through a discontinuous first-order transition 
(Fig. In the opposite case of weak anisotropy, 

K <^ J, the exchange coupling plays the dominating 
role. Accordingly the overturn of the endmost moment is 
spread over the entire stack and creates a spatially inho- 
mogeneous spin configuration (Fig. ISJiSiS In Ref. [s^ 
this mode was called inhomogeneous spin-flop phase. In 
increasing fields, a curious evolution takes place within 
this spin-flop phase, where some moments rotate against 
the applied field and change their sense of rotation at 
higher fields,--. A continuous spin- flip into the saturated 
state occurs at an "exchange" field H-^, which depends 
on the number of layers and on the anisotropy K. In the 
region of moderate anisotropy spatially inhomogeneous 
asymmetric states exist as transitional phases between 
the inhomogeneous SF and FM phases. These asym- 
metric phases arise by canting transitions, i.e., elastic 
distortions of the collinear FM phases when K/ J ^ 1. 
These asymmetric canted (C) phases can be considered 
as superpositions of ferrimagnetic states and the inhomo- 
geneous spin-fiop state (Fig. 13). This means that in these 
low symmetry intermediate C phases all the symmetries 
are broken that are broken in the corresponding SF and 
FM phases. Thus, magnetic states arising in Mills model 
comprise antiferromagnetic, spin-flop, and ferromagnetic 
phases, which exist in bulk antiferromagnets, and addi- 
tional ferrimagnetic and canted configurations. The lat- 
ter phases are imposed by the exchange cut. They are 
specific to finite antiferromagnetic layer systems. In par- 
ticular, the detailed solutions for larger N show series 
of different canted phases. The corresponding phase di- 
agrams include a large number of critical points and a 
ta ngled net of transition and lability lines (see examples 
in |47l53l54j '). The cut exchange bonds under ly this com- 
plexity as the general physical mechanism for the forma- 
tion of the various inhomogeneous magnetic states and 
their transitions into the simple collinear states in the 
limiting regions of the phase diagram for and low high 
anisotropy and for large fields. Therefore, all these phase 
diagrams have a general topology represented in Fig.|31 
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FIG. 3: (Color online) Sketch of the phase diagram for A'^ = 6 
introduces the main types of the solutions for Mills model 
Eq. J^J with arbitrary even N (A). Due to the cut exchange 
bonds the flopped states are spatially inhomogeneous and can 
exist as symmetric inhomogeneous spin-flop phase (SF) for 
low anisotropy, or as asymmetric canted phases (C) for mod- 
erate anisotropies. In strongly anisotropic systems {K > J) 
spin configurations with fiipped spins exist These collinear 
ferrimagnetic (FM) phases exist at intermediate fields be- 
tween the antiferromagnetic and ferromagnetic states. The 
spin configurations in panel (B) show for N — S the degener- 
acy of the inhomogeneous SF state for A = 8 (a) and (b); for 
A = 6 the evolution of the inhomogeneous SF in increasing 
field (c) to (e) and the two different types of FM states, where 
FMi is degenerate (f), (h) and FM2 (h) is a collinear version 
of the inhomogeneous SF state (c)-(e) in an even-odd system. 



III. PHASE DIAGRAM OF THE SOLUTIONS 

In this section we analyse the solutions for Mills model 
0, derive the regions of their existence, and conditions 
of the transitions between different phases. 



A. Spin-flop transition and solutions for 
inhomogeneous spin-flop phases 

We start with low-anisotropy systems (J ^ K). 
Here, we consider generalized models that keep mirror 
symmetry about the multilayer center with parameters 
Ji = Jn-i, Ki = Kj^+i-i, etc., for i = 1 . . . A^ — 2 or 
A^ — 1, respectively. Then, the equations, that mini- 
mize the energy of these systems, have solutions for an 
inhomogeneous spin-flop phase with the property 9i — 
— 6'Ar_i+iiSiS These magnetic configurations have differ- 
ent properties when N — 41, i.e., A^ is divisible by four 
called here even-even systems, or when N = 41 -\- 2 for 
even-odd systems {I = 1,2...). 

In low fields, Ji 3> H, the spins in the flopped state 
deviate only slightly from the the direction perpendicular 
to the easy axis 

S2j-i = 7r/2 - a2j-i, 02j = -7r/2 + aN-2j+i, (10) 
|a2j-i|«l j = l,...Af/2. 

The expansion of energy with respect to the small 
parameters ct2j-i allows to derive analytical solutions for 
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and in the even-odd case N = 10 



H/J 



^ ' KIJ=0.5 




HJ 



FIG. 4: (Color online) Evolution of the equilibrium mag- 
netic configurations 6i in Mills model for antiferromagnetic 
N = 6 layers with finite uniaxial anisotropy, K = 0.01 J (a), 
K = 0.5 J (b), in magnetic fields applied in direction of the 
easy axis. For the low anisotropy case only AF, SF, and F 
states exist (a). The points indicate the spin-fiop field -fftr 
and the spin flip field Hp. Characteristic fields are Hk, H±, 
and H* as marked in the plot. The spin configurations corre- 
spond to H = H± = 1.104 J where the magnetization in the 
2th and 5th layer is perpendicular to the applied field. For 
intermediate anisotropy (b) canted asymmetric states exist in 
the field interval Hci < H < Hc2- The spin configuration 
in (b) corresponds to the canted phase at the transition field 
Hci = 0.906 J. 



the flopped states that can be formulated for arbitrary 
models. As an illustration we write the parameters a2j^i 
for Mills model i.e., with equal constants in the even- 
even case N — 8 



2H f 11 K\ H f 6K\ , , 

2HK H f 9K 

7' "^ = 7^ + — 



_5H 
~ 27 
IH 

"^=27 
3H 
"^=27 



17 K 
2"7 
13 X 
2"7 
25 K 
2"7 



1 + — — 



1 + — — 



3H 

'"^=27 
1 H 
'"^^27 



59 K 
6"7 
37 K 



1 + — — 



(12) 



The Eqs. l(TU|l . ((TT|l . and ((T^ display the generic struc- 
ture of these solutions which apply also for generalized 
models with arbitrary parameters. The deviations of 
the magnetization direction from the directions ±7r/2 
are small of the order H/J. The corrections due to 
the anisotropy are of the order [HK)/,]"^. The devia- 
tions of the magnetization direction in the different layers 
gradually increase towards the endmost layers i = 1 or 
i = N, respectively. The dominating exchange interac- 
tions favour antiparallel ordering in the adjacent layers. 
In the flopped configurations (|11|) . (|12|l pairs of spins es- 
sentially remain antiparallel. E.g., for N — 8, the interior 
pairs S2, S3 and Sg, sy are almost antiparallel (Fig. ElB, 
panel (a)). This fact may be stated more precisely. Ac- 
cording to ini) |6'3-6l2| = \e7-ee\ = TT-{3HK)/J^, this 
means the slight deviations from antiparallel arrange- 
ment are due to a second order effect. The exchange 
coupling in such pairs is stronger than the Zeeman energy 
of the pair. This causes an interesting effect, a reverse 
rotation of the magnetization in a number of layers. In 
the multilayer with N = 8, the magnetizations S2, S7 un- 
dergo such a reverse rotation according to Eqs. Hll|l : for 
the case iV = 10 the corresponding magnetizations are 
S2, S4, S7, and sg, see Eqs. (|T^ . 

For even-odd systems, the projections of the magneti- 
zation vectors for the central layers (sjv/2, s^/2+1) onto 
the field direction are of the order H/ J, which is much 
larger than the corresponding projections of the order 
{HK)/J'^ in even-even systems (compare solutions for 
6*4,^5 in Eqs. ifTT)) . 6*5, (tT^ and spin configurations 
in panels (a) and (c) Inset B of Fig. |3| respectively). 
The solutions 9i{H) for Mills model of a multilayer with 
TV = 6 in Fig. ^ illustrate the general features of the 
field induced evolution of the spin-flop state. (See also 
the configurations in panels (c)-(e) in Inset B of Fig. O 
and the solutions for = 16 in Ref. ^] and = 12 in 
Ref. 01). An increasing magnetic field gradually slows 
down the reverse rotation of the spins with negative pro- 
jections onto the field direction. Finally at characteristic 
fields Hj^ the sense of rotation changes. In this point 

dOi/dH = 0. Another set of characteristic fields H^^ de- 
fines the points where the projection of changes the 
sign, i.e., {si{H''^) ■ H) = 0. In an increasing field these 

characteristic fields, and then , are reached first 
for the central layers and at higher fields for those closer 
to the boundaries. Fig.2|(a) displays angles 9i for an ex- 
ample of Mills model Eq. with low K ^ Q. Finally, 
for Mills model at a particular field H* the projection 
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of the magnetization onto the field direction is equal for 
all interior layers. This knot point is designated by . 
The value of H* is analytically given by 



H* = 2V3J i 1 



K 
4J 



K 
3J 



1/2 



(13) 




-1/2 



In particular, for K — the knot point is H* — \/QHy/4, 
which coincides with the value derived in Ref. 55] . For 
H > H* the positive projections of the magnetization 
onto the direction of the magnetic field decreases to- 
wards the center. The inhomogeneous spin-flop states 
for isotropic Mills models, if = 0, exist starting from 
zero field and have similar features as those for finite 
anisotropy.^^ 



B. Critical lines and multicritical points 

In this section we determine stability regions of the so- 
lutions, conditions of the phase transitions between them 
and construct the phase diagrams. Mills model Q has 
three independent control parameters K/J, H/J^ and 
iV. Correspondingly a set of {K/J, H/J) phase planes 
for different N provides complete information about the 
solutions for the model As illustration we present 
the {K/J, H/J) phase diagram for Mills model with 
= 6 (Fig. O, which containts all essential features 
of the generic phase diagram and some complications 
(Fig. I3b)), which are absent in the simplest case = 4 
as presented in Ref. The essential critical points for 
Mills models with A^ = 4 ... 16 are given in Table I. Then 
we proceed to discuss general features of the model with 
arbitrary N . 

To determine the conditions for the phase transitions 
between different spin configurations and the stability 
regions of these phases we may use standard procedures. 
The equality of the equilibrium energies of the competing 
phases yields the condition of the first-order transitions. 
The stability of the solutions {^i}i=i,...,Ar can be checked 
by writing the energy of the system for small arbitrary 
distortions, 9i = 6i + 59i, which yields the expansion 



N 



N\Oi) 



Ay = d^^N/d9^d9j . 



(14) 



The solutions 9i are stable, if all eigenvalue of the sym- 
metric matrix Apf — {{Aij)) are positive. In particular, 
for the AF phase within Mills model the matrix Apf has 
a band structure given by ^2j-i.2i-i = J + K + H, 
Mj, 2j = J + K - H iov j = I... N/2, and Ai^ i+i = 
Ai+i^ i — J for j = 1 ... A^ — 1. All other elements are 



equal zero. The determinant = det(Ajv) can be re- 
duced to the following form 



D 



N 



afJ 



(15) 



x{[i/2- {K"^ +UK + 2j'^)]DN^i 

N/2 

+ [H^ - {K^ + UK + 3J2)] Y,{-lyDN-2^} . 

1=3 

The obvious convention Do = 1 and Du = for A; < 
starts the recursion in Eq. (|16|l . Any determinant Djsi 
includes the determinant D2 for a two-layer system D2 as 
a multiplier. Thus, within the Mills model (0 the lability 
field of the AF phase has the same value for arbitrary 
values of jViiSi^Si^M^. It coincides with that for a hulk 
antiferromagnet 7?af — \J K (2 J + K) |(SJ). 

Note that this simple result for the stability limit of 
the AF phase holds only for Mills model because of its 
high symmetry. For general models (see Eq. (|17(l ^. i?AF 
is an involved combination of the materials parameters 
and depends on A^. For example, for the modified Mills 
model one derives 



Dt, = {H^ +piH^ + qi) DN-i 

N/2 



(16) 



where pi = 3J2 - ( J + K)^ - (J, + K,)^, p2 = Pi ~ J^ 
qi = qo - J'^{Js + Ks)'^, q2 = qo - ■J''', and go = [J^ - 
(J + K){Js + ks)Y. The determinants DN_2j {j > 2) in 
the right part of Eq. (|17ll are sub-determinants and do 
not include "surface" terms. 

In Fig. El the low-anisotropy region is shown by AH- 
AT-phase diagrams for A^=6 as representative for the gen- 
eral behaviour. Here, magnetic fields are given relative 
to Haf, i-c. AH = H — Hap. The stability limit for the 
antiferromagnetic phase is always above the lower stabil- 
ity limits of the symmetric and inhomogeneous SF-phase. 

Comparing the equilibrium energies in the AF and SF 
phase we determine the field for the first-order transi- 
tion between these two phases (Fig. 01. In the low- 
anisotropy limit this first-order transition line and the 
two lability fields for the stability limits of the AF and SF 
phase, respectively, are close to the value Hq = \/2 JH 
from Eq. ©. The difference between them defining 
the co-existence region for metastable states is of order 
H^{K/J) « iJo. 

Near the spin-flip transition from the SF to the F 
phase, the deviations of from the field directions are 
small {9i <C 1), and for model Eq. Q the stability matrix 

in Eq. lfT3|l becomes $Ar(6'i) = J2i j=i ^v^'.^j ^i^^ Aij 
having a tridiagonal band matrix form where nonzero el- 
ements occur only in the main diagonal and the first side 
diagonals. In particular, for Mills model A^ = H+K—J, 
Ai^i-i = Ai-i^i = J, and the spin-flip or exchange field 
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Anisotropy K/J Anisotropy K/J 

FIG. 5: (Color online) Phase diagram for Mills model with A''=6 (a). Detail of phase diagram for Mills model with A''=6 (b). 
Special canted states Cx (sketched configuration) are stable in the hatched area. The point y is a triple point, where FM, SF, 
and Cx phase coexist. Triangles designate tricritical points, where the SF-Cx transitions change from continuous to first-order 
with increasing anisotropy. 




0.1 0.2 0.3 0.4 0.5 
Anisotropy K/J 



FIG. 6: (Color online) Low anisotropy region of the phase 
diagram for Mills model with A'^ = 6. The magnetic field is 
given relative to the stability limit of the AF phase by AH = 
H — Haf{K). The canted phase Ci is stable in the darker 
(blue) area. It is metastable in the brighter (turquoise) areas. 
Thick black lines give the first-order spin-flop transition from 
AF to SF and AF to C, respectively. For anisotropies below 
point a only a first-order spin-flop from AF to the symmetric 
SF-phase exists. At higher anisotropies above the point b, 
the first-order transition is from SF to the asymmetric canted 
C-phase. For Ka < K < Kb, the canted phase exists only as 
a metastable state. For K > Kc above the tricritical point c, 
the transition from the asymmetric phase Ci to the spin-flop 
phase is flrst-order. Line a-a-b is the lower stability limit of 
the AF-phase, along the line a-b-c the continuous transition 
between C and SF takes place. Lines starting at point a are 
the metastability limits of the canted phase C. 
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FIG. 7: Hysteresis for Mills model A'' = 6 and an anisotropy 
value, where a canting instability of the spin-flop phase oc- 
curs in the metastability region Ka < K < Ki,. A resulting 
magnetic configuration of the canted Ci-state is shown. Field 
Hti for the first-order transition between antiferromagnetic 
and spin-flop-phase, and upper stability limit Haf of the an- 
tiferromagnetic phase are indicated. 

is a linear function of K 

iJ^^^ (K) = 2J + Kf-K, i/^^^ (0) ^2J + Kf,{17) 

where Kf is defined as the value of uniaxial anisotropy 
in the point / where, depending on iV, the hue -ff^^'' (K) 
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FIG. 8: (Color online) Equilibrium solutions for Mills model 
with N=16 and intermediate anisotropy K/J = 0.5. In in- 
creasing field a series of transition lead from antiferromagnetic 
to ferromagnetic phase via several canted (Ci) and (reentrant) 
spin-flop states. Left panel: magnetization of the equilibrium 
states(first-order transitions are marked by arrows). Right 
panel: spin-configurations in the canted in spin-fiop states. 



intersects line H^m = 2 J. The line H^m is the transition 
line between the ferromagnetic and ferrimagnetic phases. 
The value Kf can be derived analytically as solutions of 
the equation DN{Aij) ~ with An = Kf — J, Ai^i-i = 

= J ( Table I); H^^\o) = 2J + Kf is the spin-flip 
field for zero anisotropy (in point g, Fig.\^. In Ref. [53 
the spin-flip flelds have been calculated for generalized 
isotropic models including biquadratic exchange. 

For systems with larger anisotropics, an asymmetric 



canted phase Ci occurs flrst as a metastable state for 
Ka < K < Kf,, which can be reached by a continu- 
ous canting of the spin-flop phase. A corresponding hys- 
teresis around the first-order transition between AF and 
SF-phase is shown in Fig. lO with an example of the 
magnetic configuration Ci. This canted phase Ci is de- 
rived from elastically distorting the coUinear ferrimag- 
netic state FMi with a ferromagnetically aligned pair at 
the surface. For even larger anisotropy Kt < K, the 
Ci state becomes a stable phase of the system, which 
is reached from the AF-state through a first-order tran- 
sition. In the H vs. K-phase diagrams, the magnetic 
fields for the upper and lower stability limit of the canted 
phase Ci meet at the critical point a at {Ka, Ha) This 
point also delimits the line for the canting instability of 
the metastable SF-state. The critical line for the cant- 
ing instability of the stable SF-state ends in the criti- 
cal point {Kh, Hh) on the line of first-order transitions 
Hti- between either the AF-phase and the SF-phase be- 
low Kh, or the AF-phase and the Ci -phase above Kt- 
This point b located at {Kb, Hh) designates the lower 
anisotropy limit, where the phase Ci and any asymmet- 
ric canted phase is stable for Mills models. The topology 
of the phase diagrams in Fig. for the corner of low 
anisotropy K < Ka and fields describes the general be- 
haviour for arbitrary N. From our previous analysis, we 
have seen that no canted asymmetric phase may occur at 
low anisotropics. The first canting instability at higher 
anisotropy will occur into a phase similar to the Ci phase 
with a fiopped configuration at the surface. For Mills 
models with various TV we have numerically determined 
the low-anisotropy parts of the phase diagrams and veri- 
fied this general topology. Coordinates of the two critical 
points a and b for the canting instabilities are collected 
in Table m from numerical investigations of Mills models 
with iV = 4,6, .., 16. 

Magnetization curves corresponding to the equilibrium 
states, where the canted state Ci is a stable phase are 
shown in Fig. |7| For anisotropy K > Ki, further tran- 
sitions and critical points occur depending on N. E.g., 
the transition between the Ci-phase and the spin-fiop for 
N = 6 becomes first-order above a tricritical point c. 



C. Metamagnetism of strongly anisotropic systems 

At high enough uniaxial anisotropy only coUinear 
phases are stable. For the infinite anisotropy limit, one 
can describe the model as an antiferromagnetic chain of 
classical Ising-spins. In this limit, all coUinear states co- 
exist and transitions between them are first-order. The 
equilibrium states and their transitions are found from 
the comparison of their different Zceman and exchange 
energy. In Mills model Q only two first order transi- 
tions take place. For increasing fields these are a tran- 
sition from the antiferromagnetic (AF) state to a set of 
degenerate ferrimagnetic (FM) phases at H = J, and the 
transitions from these FM-phase into the saturated phase 
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TABLE I: List of main critical points in the field-anisotropy phase diagrams for Mills model antiferromagnetic multilayers 
with A'' = 2, 4, . . . , 16. Field and anisotropy values are given in units of J. 



N 


(K H \ 


(Kh Hh) 




^ »- e 




4 


(0.160, 0.522) 


(0.300, 0.730) 


(0.622102, 1.57956) 


0.847759 


21/2 


6 


(0.090, 0.408) 


(0.206, 0.620) 


(0.637223, 1.51922) 


0.842236 


31/2 


8 


(0.051, 0.312) 


(0.120, 0.481) 


(0.639260, 1.50798) 


0.842001 


V2 + 21/2 


10 


(0.034, 0.256) 


(0.080, 0.394) 


(0.639621, 1.50545) 


0.8419914 


V10 + 20V2/2 


12 
14 


(0.024, 0.217) 
(0.019, 0.193) 


(0.056, 0.332) 
(0.042, 0.286) 


(0.639689, 1.50486) 
(0.639702, 1.50472) 


0.841990990 
0.8419909729 


(gl/2_^2l/2)/2 

1.9498* 


16 


(0.014, 0.169) 


(0.032, 0.251) 


(0.639705, 1.50469) 


0.8419909721 


\/2 + V2 + 21/2 



"calculated from analytic expression with arbitrary precision. 
'Numerical value given instead of a long analytic expression. 



(F) &iH ^ 2 J. 

Due to the high symmetry of Mills model, it displays 
a remarkable degeneracy of the FM phase. This degen- 
eracy has important consequences for the structure of 
the phase diagram at finite anisotropy. Let us denote 
a ferromagnetic pair with configuration ]]. The two 
different antiferromagnetic domains are (AF1)=|J, and 
with reversed spins (AF2)=|t. The ferrimagnetic con- 
figuration with a flipped spin at the edge can be writ- 
ten (AFl)^/^~-'^(F), where exponents denote the number 
of repetitions for a pair. It is easy to see that all con- 
figurations of type FM„ = (AFl)^/2-n-i(-p)(^p2)" with 
n = 1, . . . , N/2 have the same energy for Mills model 
(Fig. O (b) panel (f)-(h)). There are no further ferri- 
magnetic equilibrium phases for this model. For gener- 
alized models with differing magnetic properties of indi- 
vidual layers the degeneracy of the FM phases will be 
lifted. Then, the two transitions for Mills model in the 
limit of infinite anisotropy will be replaced by sequences 
of metamagnetic transitions between various asymmetric 
collinear states. The exact sequence will be subject to 
the set of materials parameters for the individual layers. 

Towards finite anisotropy, the collinear phases will un- 
dergo characteristic instabilities were the competition of 
Zeeman energy, exchange and anisotropy will cause elas- 
tic distortions of these configurations. The stability lim- 
its for these collinear phases can be calculated from the 
analytic expressions for zero eigenvalues of their stability 
matrices A in Eq. (|15|) . These analytic expressions are 
derived similarly to the expressions for the upper stabil- 
ity limit _ffAF of the AF phase, Eqs. ijHJl, and for 
the lower stability limit or exchange field h'^^ of the 
F phase, Eq. lfT7|l . In principle, they can be evaluated 
with arbitrary precision. But, the stability limits of the 
different FM phases depend not only on N but also on 
the particular realization FM;, i.e., on the location of the 
F-pair in the chain. For the ferrimagnetic phase of Mills 
model with = 4 there is only one lability line iJpM [K] . 
It can be written in the parametric form 



FM 



J{t-K{t)), K 



Jt{t^ - 2^ + 2) 
2(^3 - 2t - 1) 



(18) 



The line Hpm{K) H18() consists of two branches meeting 
in the point {K*, H*) where K* = K{t*) = 0.62210, 
H* ^t* - K{t*) = 1.57956, {t*f = 2 + 2^/^ + 4^/3. 

For arbitrary N , the stability limits for the various en- 
ergetically degenerate ferrimagnetic phases FM; are dif- 
ferent. Here, I = 1 . . . [{N — 2)/2] and square brackets 
[x] denote the largest integer I < x. However, they dis- 
play a certain systematics for Mills models. This can 
be understood from the weakened exchange stiffness at 
the surfaces, i.e. the surface cut, which distinguishes the 
state FMi with an F-pair at the surface from all other 
realizations F; with / > 1 (see Fig. EJb) for the sim- 
plest case N = Q). The generic behaviour of these lines 
is demonstrated for the case A^ = 16 in Fig. O The 
lower branch for FMi occurs always at higher fields and 
anisotropics. This is the expected consequence of the 
surface cut. The FMi structure is more easily distorted 
and plays a special role in the intermediate anisotropy 
range. In all phase diagrams, at a certain section of this 
line a continuous transition between the FMi phase and 
the corresponding Ci-phase occurs (see Fig.l^^b)). This 
transition line at low fields and high anisotropies, starts 
at the critical end-point (ife, Hf.) on the first-order line 
between the AF and FM-phase. In the phase diagrams 
for A^ — 6, this section ends at a multicritical point x. 
For general A^, the other end of the line of continuous 
transitions between the FMi and the Ci phase depends 
on N because at higher fields other canted phases C; 
with / > 1 may intervene. Thus, the series of co-existing 
collinear states FM; at high anisotropy gives rise to se- 
ries of corresponding canted phases C; by elastic distor- 
tions. These canted phases, however, are stable in dif- 
ferent regions of the phase diagram towards intermediate 
anisotropies. Thus, the series of canted phases in increas- 
ing fields starts with Ci and, for the highly symmetric 
Mills model, it follows the sequence / = 1, . . . (Fig.|Sfb)). 

For the first canting transition between the ferrimag- 
netic state FMi and the corresponding canted phase Ci, 
along the line for the stability limit of the FMi phase, 
there is a minimum value of anisotropy Kp. The corre- 
sponding critical point (3 coincides with the upper limit of 
fields Hp for the metastability limit of the canted phase 
Ci. The values of these two characteristic points are 
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FIG. 9: (Color online) High-anisotropy phase diagram for 
Mills model and A'^=16. The degenerate ferrimagnetic phases 
are stable in the shaded area. The coUinear phases FMi and 
FM2 to FM4 are hmited by the lability hnes (1) and (2-4). 
The latter are nearly degenerate and not resolved. 

listed in Table |T| for Mills models N =4 to 16. There is 
only a weak dependence on N for the coordinates for the 
critical points (3 and e. This means, that there is also 
only a weak shift of the region of stable and metastable 
FM-states. Further, the minimum anisotropy values for 
the stability region of the other coUinear phases FM/ 
with / > 1 are always larger than K13. Thus, canted 
phases may occur in the region rather well circumscribed 
by the area shown for the simplest phase diagrams for 
= 4 and 6. Interestingly, all the stability lines for 
these coUinear phases FM; / — 1, . . . cross at the point 
X at {K^, H^) = (4/3J, 2/3J) for Mills model and arbi- 
trary N. One can show this by a similar recursion for the 
eigenvalues of their tridiagonal stability matrices as used 
by Dantas et al. to calculate Hp^p^ This point is also 
visited by the line -/Jaf {K) . For all systems with an even- 
odd number of layers iV=6, 10, etc., the point a; is a spe- 
cial multicritical point, where the coUinear analogue of 
the symmetric spin- flop phase FM[(jv-2)/2] becomes de- 
generate with the inhomogeneous symmetric SF- state. 
Thus, the first-order transition line between the corre- 
sponding canted phase C[(Ar_2)/2] and the SF-phase ends 
in X as in the phase diagram for = 6 shown in Fig. 
This further degeneracy of Mills model yields some sta- 
bility and simplicity of the general features of its phase 
diagrams. 

As a caveat, we finally have to mention additional low- 
symmetry phases which cannot be foreseen from the con- 
siderations on the stable coUinear phases AF, FM;, and F 
and their clastic distortions in the phase diagram. Such 
an intermediate phase Cx appears already in the phase 
diagram of Mills model for = 6 in the region between 
the FM; phases and the region of the stable SF-phase. 
The Gx phase can be derived from an elastic distortion of 
the coUinear phases TiTTTTj which do not arise as stable 
states. Similar low-symmetry canted phases also appear 



in phase diagrams of Mills model for iV > 6. For gen- 
eralized models, the region of stability of these canted 
low-symmetry phases will strongly depend on details of 
the materials parameters of type Eq. Q . Further energy 
terms and/or disorder in magnetic parameters of the lay- 
ers can stabilize further phases in the intermediate region 
of the phase diagram. In the high-anisotropy limit of such 
generalized models, cascades of metamagnetic transitions 
between ferrimagnetic coUinear phases exist. From these 
states various canted phases can be derived in the range 
of intermediate anisotropics Kij Ji ^ 1. The competition 
between all these phases will lead to very complicated 
magnetic phase diagram 



IV. REORIENTATION IN WEAKLY 
ANISOTROPIC MULTILAYERS 

A. Surface and volume interactions 

In this section, we re-analyze the magnetic states of 
an antiferromagnetic multilayer stack, and their evo- 
lution in an applied field, by analytical methods. In 
the limit of weak anisotropy, the micromagnetic en- 
ergy JQ) can be represented by a system of interacting 
dimers and by a continuum form. Thus, we consider 
the model for weak anisotropics, Ki, K[ ^ Ji. Fur- 
ther, we only study systems with coUinear antiferromag- 
netic (zero-field) ground state. Thus, the strengths of the 
biquadratic exchange constants is limited to the range, 
< Ji < Ji/2.^^ First, we rewrite the general energy 
Eq. (^) in this limit. The resulting expression allows to 
recognize the main effects expected in this limit without 
explicit calculations. We group the moments of the su- 
perlattice into pairs as in a two-sublattice antiferromag- 
netic system. Starting from the first layer we combine the 
N moments along the antiferromagnetic chain into N/2 
pairs (s2j_i,S2j) with j — \,..,N/2. For each of these 
pairs, we introduce the vectors of net magnetization Mj 
and the staggered vector Lj 

S2,_i=M,+L„ S2j=M, -L,. (19) 

These transformations are similar to those applied for the 
two-layer system, Eqs. Q and {Tj). From |si| = 1 follows 
that Mjjr L"] = 1 and Mj • = 0. The energy 9^ 
of Eq. (^3) can be rewritten as a function of Mj and the 
angles <j)j between n and unity vectors \j — Lj / 1 Lj | and 
expanded with respect to the small parameters Mj <C 1. 
Omitting constant and higher order terms, one derives 

N/2 N/2 N/2 

Qn = ^AjM2-2^(H-Mj)-^^iCos2(/)j 

j=i j=i j=i 

N/2-1 

+ 2 E J2,(0,+i-0,f + S(M„0,), (20) 
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wherein, = {K2,-i+K2j)/2^K'^^^^-{K'^^_^+K^^^)/2, 

= 2(J2j_l-2J2i-l) + (J2j-2-2J2j-2) + (J2j-2J2j) for 

j^- 2, . . . , N/2 - 1, and = {K^ + K^)I2 -K[- K'J2, 
Kn/2 = {Kn-i+Kn)/2-K'j^_^^K'j^_^/2; Ai = 2(Ji- 

2Ji) + (J2 — 2J2), A^/2 — 2(Jn-1 — 2J]y^i) + {Jn-2 — 

2Jn-2)- Finally, the last expression in Eq. 1)20(1 collects 
terms that are linear with respect to ((/>j+i — 4>j), 

S(M„0,) = (21) 
jv/2-1 

- (-^2^ - J2j ) (Mj + Af,+i ) - 0, ) . 

An independe nt min imization with respect to Mj (see 
details in Ref. |25l66l |'l yields 

M, =A7i[H-(H.l,)y . (22) 

It follows directly from if^ that — HAj^ sin(0j — V); 
where ip is the angle between the field H and the easy 
axis n in the model of Eq. . 

The independent minimization with respect to Mj is 
justified because the exchange interactions dominate the 
energy and pairs of neighbouring moments do not devi- 
ate strongly from antiparallel orientation in the limit of 
weak anisotropy and fields. This establishes the relations 
Eq. H22|l between the components of the net magnetiza- 
tion and the orientation of the staggered vector. In other 
words Eq. 1)22(1 fixedly connects the net magnetizations, 
as auxiliary degrees of freedom, to the vectors Ij and the 
applied field. This approach reduces the chain of N mag- 
netic moments s,; into an equivalent system of A'^/2 unity 
vectors Ij. Each site of this chain corresponds to a two- 
sublattice antiferromagnet or a dimer. Substituting l(22() 
into Eq. ((20(1 we obtain the following expression for the 
energy of these interacting dimers 

N/2 

Qn = -^^j cos 2((/)j - ^j) (23) 

^ 7V/2-1 

+ 2 E ■^2,(<^, + l-0jf+S(0j), 
J = l 

$j = ^ a/ (i72 cos 2^j-H^f + i?4 sin^ 2^/^ , (24) 
2Aj 

tan2<^j = sm2ilj / {H^ cos2ij - Hf), (25) 
Hj = \/^, (26) 

where 

N/2-1 

E{<t>j) = -H f^,sin(0,~7A) (</.,+! -0,), (27) 
i=i 

^3 = J2d^]'+^]li) ■ (28) 



The minimization with respect to Mj according to 
Eq. ((2211 and the representation of the energy by the 
form ((23(1 generalizes simplified dimeriza tion tr ansforma- 
tion that have been considered in Refs. '54'67l. 

The energy of the interacting dimers Eq. ((23(1 includes 
first the sum of their "self" -energies, then an exchange 
stiffness energy given by the term quadratic with respect 
to differences {4>j+i — 4>j), and a specific energy contribu- 
tion S((/)j), defined in Eq. ((?7jl . The terms in S(<^j) arise 
due to the variation of the magnetic parameters along 
the chain. The energy S(</)-,) can be written in the form 
of a "Zeeman energy" for the staggered magnetization 
vectors 

S(0,) = -[r!i(H.li)-f}^,/2_i(H.l^./2)] (29) 
Ar/2-2 

- (%+i-%)(H-l,+i) . 

The dimensionless coefficients VI j are ratios of exchange 
constants defined in Eq. ((28(1 . The first two terms in 
Eq. I(29() involve the endmost dimers, i.e., they have the 
character of a surface energy, which is imposed by the 
exchange cut. The sum in ((29(1 describes similar "inter- 
nal" contributions arising due to any variation of the ex- 
change couplings along the antiferromagnetic chain. This 
energy contribution disappears in models with equal ex- 
change interactions in internal layers as in the regular 
and modified Mills models Eqs.|3Jl and (O, respectively. 

B. Physical mechanism of the "surface spin-flop" 
phenomena 

The magnetic energy of the low-anisotropy antiferro- 
magnetic multilayers in the form of interacting dimers 
((23(1 elucidates the competing forces responsible for the 
field-driven reorientation processes. Let us compare en- 
ergy l(23() with that of an isolated dimer. For a localized 
pair of S2j-i and S2j spins (i.e. j-th dimer) a minimiza- 
tion via Eq. (|^ yields 

e^^' = -<i>f cos2(0,-0f ), (30) 
with an amplitude factor 

= ^ {m cos 2i; - Hj) ' + iJ4 sin^ 2^, ,(31) 

and a "phase" 

tan 20^°^ ==ij2sin2V'/(i?^cos2V'-ff|), (32) 

where = [(J2j_i-2J2,_i)(i^2,_i+i^2,+2i^^,_i)]i/2. 
The energy in the form 1(30(1 coincides with that of a two- 
sublattice antiferromagnet and is a generalization of the 
model Eq. Q . The Eq. (j32Jl for the phase 0^"^ is known 
as Neel equationP'^ It determines the equilibrium states 
of the antiferromagnet </> = +TTk {k — 0,1, . . .). The 
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amplitude from Eq. 1)31(1 equals the potential bar- 
rier between the wells at (j) = (f>i + irk. A magnetic field 
applied along the easy axis reduces the potential barrier. 
When the field reaches the threshold field Hj it causes 
the spin-flop transition. For dimers incorporated into 
the interacting chain the parameters of the self-energies 
are modified due to the exchange coupling and additional 
anisotropy contributions as seen by comparing Eqs. H23|l - 
(1231 with Eqs. ((SOIl-lEll- Therefore, within the chain 
the dimers have different threshold fields and they are 
elastically coupled with neighbouring pairs. Due to the 
couplings the flopping of the individual dimers in their 
individual threshold fields Hj Eq. H23I) are hampered. In- 
stead the chain only can transform into the flopped phase 
when the flopped conflgurations are energetically advan- 
tageous throughout the whole system. Generally the dif- 
ferences of the dimer self-energies along the chain causes 
spatial modulations of any noncoUinear magnetic states. 
The inhomogeneous spin-flop and canted phases in Mills 
model (Fig. O exemplify such spin-configurations. 

However, the energy contribution S((/)j) Eq. H21|l pro- 
vides another mechanism of magnetic-field-induced reori- 
entation imposed by the variation of the exchange inter- 
actions along the antiferromagnetic superlattice, in par- 
ticular by the exchange cut at its ends. This mechanism 
is due to the influence of the linear energy terms (|27|l . 
which favour a rotation of the staggered vector. As can 
be seen from the equivalent Eq. ((29|l , an instability of the 
coUinear configuration is caused by the "Zeeman terms" 
that are linear in the staggered vectors Ij . Generally, the 
first term related to the two surfaces will dominate, and 
this difference will favour a rotation of these Ij. This 
enforces an inhomogeneous spin-flop phase above a cer- 
tain field. As was shown in the previous section in strong 
anisotropy systems the exchange cut leads to flips of the 
magnetization and a transition into coUinear FM phases, 
which are also inhomogeneous. In low-anisotropy sys- 
tems, under the dominating influence of the exchange 
interactions, the influence of this "local" defect is spread 
along the chain and stabilizes a spatially inhomogeneous 
structure. 

Thus, we have the following important conclusion. 
There are two different mechanisms of the field-induced 
reorientation in finite antiferromagnetic superlattices: (i) 
One of them is connected with a switching of the poten- 
tial wells in the energy of the uniaxial antiferromagnetic 
units. This mechanims is similar to the usual field-driven 
spin reorientation in (low-anisotropy) bulk antiferromag- 
nets and two-layer systems. Therefore, it is a common 
spin-flop mechanism, (ii) The other mechanism is due 
to variation of the exchange coupling along the superlat- 
tice and, in particular, the exchange cut at the end of 
the stack. This type of reorientation transition can only 
exist in finite antiferromagnetic superlattices and has no 
analogue in bulk antiferromagnetism. 

The interplay of these two mechanisms rules the forma- 
tion and evolution of the magnetic states in the system. 
Depending on the values of the material parameters dif- 



ferent types of magnetization processes can be realized 
in the general model In the low-anisotropy systems, 
owing to the dominance of exchange interactions, it is 
the second effect due to the cut exchange at the surfaces 
that dominates the field-driven reorientation transition. 

As important cases for applications, we consider in 
more detail the highly symmetric Mills models Eqs. Q 
and (|SJ|. Both models are composed of identical in- 
ternal layers. For the modified Mills model Eq. (O 
Ai = Ajv/2 = 2{2Js + J)/{1 + Cs), and A^ = 4J for 
j = 2, . . . , N/2 - 1. For the regular Mills model Eq. @ 
Aj = 4J for j = 1, . . . , N/2. The energy ^ with ip = 
reduces to 



TN 



E ^ oj ^^ cos2,/>, 

J=2 



(33) 



J 



N/2-1 



J2 (</'j+i-'^j)' + $«, 



where = \/ AJK. The last contribution in ((33() is due 
to the finite length of the chain. It involves only the two 
last dimers at both ends of the chain, i.e., it represents 
the specific surface effects for the finite antiferromagnetic 
stack. For the regular Mills model the isolated dimer 
energy l(5n|l reduces to the form of Eq. ||HJl, and the surface 
energy becomes 



7 
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1 



H'^ - Hi) (cos 201 -I- cos 2(t> 



N/2) 



H 



s 



H (cos 01 — cos cj)N/2) 

H 

V3JK. 



+ (COS02 - COS0Ar/2_i) , (34) 

(35) 



In the case of the modified Mills model, the contribu- 
tion <I>, has the modified form 



- \ (cos 201 -I- cos 20^,/2) 



8(2 J, + J) 
—UiH (cos 01 — cos 0^/2) 

-t-t/2i? (COS02 - COS0Ar/2-l) ■ 



(36) 



The threshold field Hs and the coefficients, Ui and [/2, 
are defined by 



U2 



^2{Ks + K){2J, + J) 

(1 + G) 
J(6C. - l)-H2J,(4Cs -3) 

4(2 J, + J) 
J(2C. + 1) - 2J, 
4(2 J, + J) 



(37) 
(38) 
(39) 



The threshold field for the endmost dimers is Hs = 
■\/3/ 2Ho and for internal dimers Hb = \/2Hq. Because 
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the spins in the chains have additional exchange cou- 
pHngs, these thresholds are larger than the threshold field 
Hq for an isolated pair. This reinforcing by exchange 
stiffness for the bound dimers increases the values of the 
threshold fields for the coupled chain. 



C. Multilayers with large N. Continuum model 

The limit of multilayers with large N and the limit of 
infinite antiferromagnetic systems is best discussed go- 
ing over to a continuum description. For the regular 
Mills model ^ with arbitrary N the transition into the 
flopped state occurs closely to i?tr ~ ^^o given by Eq. 0, 
sufhciently below the dimer threshold fields Hs and Hb 
(Fig. Ill|l . This means that this transition is imposed 
by the exchange cut. These results can be easily under- 
stood from the solutions for the spin-configurations in 
the flopped state Eqs. (HUJl, HH), and The magne- 

tization vectors for all internal layers can be combined 
into pairs with antiparallel magnetizations. The system 
effectively behaves as an isolated dimer consisting only of 
the endmost spins with energy Q). Correspondingly the 
flopping field equals the threshold field Hq of this isolated 
pair. This result is common for systems with arbitrary 
values of N. 

However, above the critical field Hq the evolution of the 
system remarkably changes with increasing N (Fis. Il()|l. 
The multilayer with low N are characterized by a large 
magnetization jump at the transition field and a nearly 
linear increase of the magnetization up to the flip field. 
With increasing N the magnetization jump at Hq gradu- 
ally decreases. Concurrently, a steep section of the mag- 
netization curve is found around fields with values close 
to Hb (see Eq. Finally, for ^ 1 the magnetiza- 

tion curves develop a strong kink around this value. 

The magnetic-field driven transformation of the dimer 
energies in Eq. (|33|l explains this peculiar behaviour. 
In the transition field Hn ~ Ho the dimer self-energy 
terms in (I33f) still favour the antiferromagnetic mode 
(lj ||n). The threshold fields are exceeded at higher fields 
{H > Hs for endmost and H > Hb for internal dimers). 
In superlattices with few layers the endmost and internal 
dimers give comparable contributions to the magnetic en- 
ergy. The difference in their internal energies suppresses 
drastic reorientation effects at the threshold fields Hs 
and Hb- With increasing number of layers the relative 
energy contribution of the internal dimers for the total 
energy (|20|l gradually increases. Then, the magnetic en- 
ergy of the internal layers plays the dominant role in the 
formation of the flopped configurations. Thus, the end- 
most dimers docs not hamper the reorientation effects 
in the threshold field Hb- Below the threshold field, 
H < Hb, the antiferromagnetic phase with l||n corre- 
sponds to the minima of the internal dimers and the 
inhomogeneous spin-flop phase consists of two antifer- 
romagnetic domains with antiparallel orientation of the 
staggered vectors fFig. Ill|) . These two regions are sep- 
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FIG. 10: (Color online) Magnetization curves for Mills model 
with low anisotropy K/ J— 0-01 and large N in the vicinity of 
the "bulk spin-flop" . 



arated by a 180-degree domain wall with flopped spin 
configurations in the center (1-Ln) 

For H > Hb the potential wells for the internal dimers 
switch into (1-Ln) configuration. Around the field Hb the 
center of the domain wall gradually extends and sweeps 
out the regions with antiferromagnetic configuration to- 
wards the surfaces of the stack. This drastic transforma- 
tion between the two configurations within most of the 
bulk of the antiferromagnetic stack causes a prominent 
anomaly of the magnetization curves near the field Hb 
fFig. I10|l . Above Hb the net magnetization Mj develops 
two symmetric maxima close to the surfaces, which may 
be observable in experiment. 

Asymptotically with N oo, the magnetization curve 
approaches that of the usual spin-flop in a bulk uniaxial 
antiferromagnet. But, this reorientation occurs within 
the same magnetic phase, and no phase transition is con- 
nected with this process. Rather for any finite value of 
N the phase transition still occurs between the antiferro- 
magnetic and inhomogeneous spin-flop phase in the crit- 
ical field Htr as a first-order process. As was mentioned 
above for large N this field-induced phase has the char- 
acter of a domain wall between two antiferromagnetic 
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states. Non-collinear states arise in the central region 
of the stack, where the small total magnetization of the 
configuration is concentrated. For fixed small anisotropy 
the transition is accompanied by a jump of the magneti- 
zation at the transition field i?tr- The magnitude of this 
jump decreases with the numbers of layers N (Fis. I10|l . 
Hence, at first glance we have a paradox phase diagram 
for Mills model: a drastic field-driven change of the mag- 
netization is not related to a phase transition, while a real 
phase transition is noticeable only by a small jump of the 
magnetization, that vanishes for large N. However, this 
has a clear physical foundation because the transition at 
Htr is related to the surface effect and its visible effects 
should vanish for — > oo, whereas the crossing-over to- 
wards the flopped configuration in the "bulk" of the mul- 
tilayer stack should approach a true spin-flop transition 
for N = oo. 

A transition into the inhomogeneous spin-flop state 
means that the free boundaries cause an inhomogene- 
ity far in the interior of the flnite system. Close to the 
boundaries the magnetic configuration resembles that of 
the two antiferromagnetic coUinear domains. This struc- 
ture is consistent with the properties of semi-infinite an- 
tiferromagnetic chains described by Micheletti et alifi In 
the phase diagrams for these systems (even in the large 
anisotropy limit) a highly degenerate phase occurs, where 
a localized inhomogeneous configuration is situated at ar- 
bitrary distance from the surface.*'^ For finite antiferro- 
magnetic chains with weak anisotropy, the mutual influ- 
ence of both surfaces will determine a unique state with 
180 degree wall-like configuration in the center. Gen- 
erally, such a symmetric configuration will be found for 
antiferromagnetic layers, when the core of this configu- 
ration is wide enough to interact with both surfaces. For 
the finite systems, the simple structure of the phase dia- 
gram, showing only a SF phase with solutions preserving 
mirror symmetry about the center of the layer in inter- 
mediate fields between the AF phase and the saturated 
F phase, is restricted to low anisotropy systems. At size- 
able anisotropy, the asymmetric canting complicates the 
phase-diagram. 

The energy of the modified Mills model Eq. Q pro- 
vides a simple way to introduce a continuum form of 
energy For TV ^ 1 and low anisotropy K 1 the 

energy H2U|I can be converted to 



AW ^ W - l^AF = 



d/2 



-rf/2 



A 
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i=±d/2 



(40) 



with the exchange constant A — J a/2. The multilayer 
thickness is d = Na with a the "periodicity" length of the 
multilayer. The zero of energy scale for AVF is shifted to 
the energy Waf of the antiferromagnetic state (0 = 0). 
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FIG. 11: (Color online) Evolution of configurations for Mills 
model with low anisotropy K / ,J=Q.G1 and A'^ = 100 in the 
vicinity of the "bulk spin-flop": (a) local net magnetic mo- 
ments Mj (b) orientation of local antiferromagnetic vectors 
(according to the Eq. 11911 . Inset shows average magnetiza- 
tion m vs. field. 



The last term is the surface energy given by 



4(2 J, + J) 



1 



sgn(0 (C.-§)^f[cos</)(^)-i; 



(41) 



Eq. (|40|l describes the energy of a plate of thickness d 
for a bulk easy-axis antiferromagnet with the spin-flop 
field H^. It is the continuum counterpart of the dis- 
cretized model Eq. The surface energy Ws H41|l 

includes a common antiferromagnetic contribution (the 
first term) and a specific Zeeman energy imposed by the 
exchange cut. Due to mirror symmetry of the inhomoge- 
neous spin-flop phase the boundary conditions arc (pi = 
(/)(^ = —d/2) and 02 = 0(C = rf/S) = vr — 0i. By solving 
the Eulcr equation for the energy functional H40|l with the 
boundary condition 2A{d(f)/ d£/)\^^^^/2 ~ — 9VFs(01)/30 
one obtains a set of parametrized profiles 0(^,0i). The 
further optimization of the energy with respect to the 
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parameter yields the equilibrium distribution of the 
staggered vector within the multilayer of finite thickness. 
These solutions can generally be written as elliptic func- 
tions. 

In the limit of infinite thickness, N oo the boundary 
values of (jj correspond to the configurations in antipar- 
allel domains of the antiferromagnetic phase, (f>i(—oo) 
=0, 02 (oo) = TT. In this case the variational problem 
for the functional (|40|l is equivalent to that of an iso- 
lated magnetic wall with "effective uniaxial anisotropy" 
K* = {H^ - 2H^)/{16A). The corresponding analyti- 
cal results for the wall parameters have been derived by 
Landau and Lifshitzfii. For Mills model this solution has 
been analysed ini4. 

This limiting case of the infinite chain provides a sim- 
ple physical explanation of the phase transition into the 
inhomogeneous spin-flop phase. According to Eqs. (|40|l . 
(|41|l the flop of the surface staggered vector (/)(oo) — 
into antiparallel position '/>(oo) ~ tt yields a gain of sur- 
face energy AWs = — (2Cs — l)H. By this process a 
domain wall is generated which requires a positive en- 
ergy contribution A^J AK* — ^ — . The balance 
between these competing energy contributions is reached 
at the critical field 

Hu = i/o/v/2Cf - 2C + 1 . (42) 

In particular, for the regular Mills model = 1 and 
the transition field in Eq. (|42|l equals Hq . This simple en- 
ergetics allows to formulate a clear thermodynamic rea- 
son for the transition into the inhomogeneous spin-flop 
phase provided by the exchange cut. The anti-aligned 
magnetization vector at the non- compensated surface is 
overturned and reduces the Zeeman energy at the expense 
of the formation of a planar defect in the center of the 
superlattice, which has the character of a domain wall. 

Because the energy gain in the inhomogeneous spin- 
flop phase, AVFs((s) = — (2Cs — l)-ff, is proportional to 
the non-compensated magnetization of the surface layer, 
this transition into the flopped state strongly depends 
on the net magnetization of the endmost layers. Partial 
compensation of the surface magnetization is tantamount 
to a reduction of the parameter C,s- This reduction de- 
creases AT/Fs(Cs), and increases the critical field i?tr (021 • 
For C,s = 1/2 the energy gain /S.Ws{C,s) equals zero and 
i?tr reaches the threshold field H^- 

We come to an important conclusion. The exchange 
cut provides the stabilization of the flopped phase only 
under condition of sufficiently strong surface magnetiza- 
tion. In the regular Mills model the magnetizations of 
the layers are assumed to have fixed values. The proper- 
ties of the endmost layers j — 1 and N are described by 
the same integral phenomcnological parameters as the 
layers j = 2, . . . , iV - 1 in the "bulk". Only the ex- 
change cut reflects the confinement of this antiferromag- 
netic system. In the continuum limit (|40(l the surface cut 
is represented by the surface contributions in the energy 
that describe the effective coupling of the staggered vec- 
tor to the applied field. The existence of surfaces with 



net non-compensated magnetization is justified for anti- 
ferromagnetic superlattices with in-plane magnetization. 
Strong intra-layer exchange interactions and weak stray 
field effects favour single domain states of the individual 
ferromagnetic layers in these multilayer stack and at the 
surfaces. In other systems, various mechanisms can cause 
reductions of the non-compensated magnetization at the 
surfaces, such as crystallographic and magnetic imperfec- 
tions, formation of antiferromagnetic multidomain states 
etc. A reduction of the total magnetization in the surface 
layers can strongly reduce the non-compensated magne- 
tization and suppresses eventually the formation of the 
fiopped states. For the continuum model H40|l this occurs 
for Cs < 1/2. 

The surface energy (|41ll also includes the first term that 
has the conventional form of a (local) antiferromagnetic 
unit with a modified threshold field Hs from Eq. I|37|l . 
According to many experimental observations, the mag- 
netic parameters Js^ Kg and can be strongly modified 
by surface-induced interactions, see, e.g., Refs. |6Cl69l |. 
Correspondingly Hs in l|41|) can vary in a broad range. 
Generally, considering models with modified magnetic 
surface properties, the volume energy (|40|) and the sur- 
face energy (|41|) may favour different magnetic configura- 
tions in certain intervals of applied magnetic field. This 
competition can stabilize inhomogeneous phases with 
continuous rotation of the magnetization vectors along 
the thickness. The occurrence of such twisted states un- 
der pinning (or anchoring) influence of the surfaces is a 
rather general effect in orientable media. In particular, 
they are known to occur in various classes of liquid crys- 
tals as the so-called Freedericksz effec1i22i2i and in ferro- 
magnetic materials. ^^'''^■^'^ Spiraling in exchange spring 
magnets and exchange bias systems also belongs to this 
class of phenomena/'^ The phenomcnological theory of 
such states in antiferromagnetic nanolayers has been de- 
veloped in Ref. It was shown that non-coUinear 
twisted phases can arise as solutions for magnetic states 
under anchoring-effects at the surfaces. In contrast to in- 
homogeneous spin-flop states stabilized by the exchange 
cut the twisted phases arise due to pinning or distortive 
effects of surface-induced interactions on the magnetic 
states. Future analysis of generalized Mills models should 
concentrate on the combined effects of these surface in- 
teractions. 



V. DISCUSSION 

A. Solving the "surface spin-flop" puzzle 

The exchange cut is the primary driving force that 
causes the specific reorientation effects in antiferromag- 
netically ordered multilayers and stabilizes the unique 
magnetic states unknown in other classes of magnetic sys- 
tems. The pioneering studies by Mills and co-worker9^ii4 
have introduced the notion of a surface-induced insta- 
bility in confined antiferromagnetai^ and of the novel 
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reorientation effects in antiferromagnetic superlattices^. 
Mills formulated the basic idea that, in a confined antifer- 
romagnet, uncompensated surface magnetization causes 
the instability of the coUinear state in the applied mag- 
netic field quite below the common (bulk) spin-flop. This 
transition from the antiferromagnetic to a "surface spin- 
flop" state should result in flopping a few layers of spins 
near the surface, i.e., they would turn by nearly 90 
degreejii This picture was improved and detailed by Kef- 
fer and Chowfi^, and supplemented by results of numer- 
ical simulationsi^*S4^ii& It constitutes the recent scenario 
of a "surface spin-flop" . According to this picture the 
flopped states are nucleated initially at the surface and 
in increasing flelds this surface state moves into the depth 
of the sample as an antiferromagnetic domain wall:^^''*^ 

"When the external field exceeds the surface spin-flop 
field, the surface moment, initially antiparallel to the 
field, rotates nearly by 18(F . In effect, a twist is applied 
to one end of the structure. A domain wall is then set 
up, in an off center position in the finite structure. [. ■ ■ ] 

With further increase in field, the wall undergoes a series 
of discontinuous jumps, as it migrates to the center of the 
structure. [. ■ . ] The domain wall becomes centered in the 
structure, and then with further increase in field broad- 
ens, to open up as a flower to evolve into a bulk spin-flop 
state. The angle between the spins and the external field 
is less at and near the surface than in the center of the 
structure. 

The detailed investigations of Mills model (j2Jl necessi- 
tate a considerable revision of the surface spin-flop sce- 
nario expanded inSii4 and some other papersSiiSii^. This 
scenario of the surface spin-flop confuses three different 
types of the solutions for Mills model (Q, see the phase 
diagrams in Figs. |3| and [S] The "flopping of a few lay- 
ers of spins near the surface" is inspired by solutions for 
ferrimagnetic and canted phases in systems with sizeable 
anisotropy. The picture of the domain wall movement 
into the center of the superlattice "in a sequence of dis- 
crete hops"— is related to cascades of phase transitions 
between canted phases in superlattices with moderate 
anisotropy (Figs. Finally, the flower- like broadening 
of the centered domain wall poetizes the evolution of the 
inhomogeneous spin-flop phases in flelds higher than Hb 
(Fig.Hni. 

Thus, the common scenario of the surface spin-flop 
combines elements that belong to different solutions for 
different values of the control parameters {K/J, H / J, N) 
of the model ©• In Ref.^^, it was shown that the equa- 
tions minimizing energy (|2Jl , as well as general models , 
do not include solutions for surface-conflncd (localized) 
states, which were assumed to occur at a "surface spin- 
flop" transition" . These models own only well-deflned 
"volume" phases and transitions between them. 

The term "surface spin-flop" designates the reorienta- 
tion into the inhomogeneous spin-flop phase at i?tr (line 
a — & in Fig.[SJ). Therefore, it is a double misnomer. This 
transition does not take place at the surface because it 
involves the reorientation of all spins along the superlat- 



tice, i.e., it has a "volume" character. And it is not a 
proper spin-flop because it is induced by the exchange 
cut rather than a switching of the potential wells as in 
the common spin-flops in bulk antifcrromagnets. 



B. Notes on the experimental observations of 
surface spin-flop phenomena 

The concept of a "surface spin-flop" is commonly ap- 
plied to analyse experimental results in antiferromagnetic 
superlattices^*i2iiL2i. However, the application of an er- 
roneous concept is dangerous. In particular, quantitative 
conclusions about magnetic materials parameters from 
the observed reorientation transitions can lead to wrong 
results. The lability field of the antiferromagnetic states 
Haf plays the prime role in the common "surface spin- 
flop" scenario. Because the surface spin-flop was believed 
to arise as a local surface instability of the coUinear phase 
exactly at the critical held Haf, this was considered as a 
transition held into the surface spin-flop stat o^'^i^^ . In re- 
ality a (volume) first-order transition between the antifer- 
romagnetic and inhomogeneous spin-flop phases occurs 
at Hti- (e.g., the line a-b in Fig. [SJ, which is lower than 
Haf (line aA) and larger than another lability field iJsF 
(line a-a). The interval HgF < H < Haf is a metasta- 
bility region of these competing phases (Fig. (SJ. In the 
low-anisotropy limit the metastablity region is extremely 
small and these characteristic field are all close to the 
value Hq from Eq. pOfl . In the limit of large anisotropy 
the lability field Haf is much larger than the transition 
field between AF and FM phases (Fig. EJ. 

The "bulk" spin- flop field is also considered as impor- 
tant element of the common scenario. Starting at Haf 
the expansion of the surface spin-flop phase is completed 
in flelds exactly equal to the value of the spin-flop transi- 
tion in a "bulk" antiferromagnet having the same values 
of the magnetic parameters as in model For low- 
anisotropy systems this field equals Hb and is times 
larger than the "surface spin-flop" Hq. This field corre- 
sponds to the threshold field for the internal dimers as 
given in Eq. (|33|l . In systems with large numbers of layers 
N, a strong reorientation of the fiopped states occurs in 
the vicinity of this field. No phase transition is connected 
with this process, however, it is marked by a noticeable 
anomaly of the magnetization curve (Fig. I1U|) . Thus, the 
magnetization curve anomalies are connected with the 
transition into the flopped state at Hti ~ Hq and with 
a smooth reorientation near Hb that does not involve a 
real transition. The ratio iJe/^^tr is about y/2. A simi- 
lar anomaly within the spin-flop state is also observed in 
systems with rather large anisotropy, where the spin-flop 
phase is preceded by one or several canted phases. How- 
ever, a glance at the phase diagram, e.g., Fig. O shows 
that there is no simple quantitative relation between the 
various reorientation anomalies observable in such mul- 
tilayer systems with sizeable anisotropy. 

Magnetic-fleld-induced reorientation transitions were 
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investigated in high-quality Fe/Cr(211) antiferromag- 
netic superlattice o^i^° . In Ref. [3 magnetization curves 
for Fe/Cr (211) superlattices with strong uniaxial 
anisotropy were measured by a SQUID magnetometer 
and by longitudinal magneto-optic Kerr effect. The 
magnetization curves for both investigated multilayers 
with even number of layers Cr(100)/[Fe(40)/Cr(ll)]22 
and Cr(100)/[Fe(20)/Cr(ll)]2o demonstrate close corre- 
spondence to theoretical results for Mills model. Accord- 
ing to^ the values of the antiferromagnetic coupling be- 
tween the layers is JMq = 0.275 erg/cm^ and uniaxial 
anisotropy is KMq — 0.06 erg/cm^, The ratio K/J = 
0.22 shows that these multilayers belong to the systems 
at intermediate anisotropy in the phase diagram that dis- 
play cascades of canted phases. Indeed, the characteristic 
anomalies in the field derivatives of the magnetization re- 
veal a series of such reorientation transitions. The asym- 
metric character of these transitions is demonstrated by 
Kerr measurements (see Fig. 3 (b) in Rcf. 5J). A cascade 
of canted phases has also been observed in another Fe/Cr 
systemia. In this paper a Cr(100)/[Fe(14)/Cr(ll)]2o sys- 
tem has been investigated with JMq = 0.405 erg/cm^ 
and KM^ = 0.06 erg/cm^. The ratio K/J = 0.15 means 
that this superlattice also evolves in the applied field via 
a cascade of canted phases. It should be noted that the 
first-order transitions between these different magnetic 
phases allows fo r phas e co-existence in rather wide field 
ranges (see Refs. lKHIHi V Ah these processes may involve 
multidomain states. In the case of multilayers with ef- 
fective in-plane anisotropy, the stabilization of domain 
structures will be subject to imperfections. In particu- 
lar, interface roughness will lead to magnetic charges or 
leaking dipolar stray fields. The corresponding domain 
structures is determined by the defect structure of the 
multilayer and will have an irregular appearance in gen- 
eral (see, e.g., chap. 5.5.7 in Ref. 0). 

In contrast, multilayers with perpendicular 
anisotropics constitute a novel class of artificial confined 
antiferromagnets, where well-defined and regular domain 
structures such as stripes or bubbles can be found. These 
are layered systems as antiferromagnetically coupled 
multilayers [CoPt]/Riii, or Fe-Au superlattices^i. These 
strongly anisotropic systems correspond to the right side 
of the phase diagrams in Figs. |21 Eland are characterized 
by a number metamagnetic jumpsi^. Due to strong de- 
magnetization fields the magnetization processes in these 
superlattices are accompanied by a complex evolution of 
multidomain states2i2iSSiS4. Artificial layered systems of 
this kind with a controlled variation of distinct magnetic 
states in vertical direction can be considered as artificial 
metamagnets^. 

VI. CONCLUSIONS 

In this work, we have provided a complete solution for 
the basic micromagnetic model of an antiferromagnetic 



superlattice with ideal non-compensated surfaces under 
a field along the easy axis. We have shown how one 
can systematically enumerate and describe the magnetic 
phases and their transitions for such structures. The puz- 
zle of the variable appearance of "surface spin-flop" phe- 
nomena has been resolved by the re-construction of the 
phase-diagrams and of the limiting cases for this model. 
To this end various methods had to be introduced that 
can be used for generalized models. Analytical tools can 
be efficiently used for all coUinear or highly symmetric 
phases, and for the case of weak anisotropics. Exten- 
sions as given by the models Eqs. ^ ^©i that include 
further magnetic coupling terms, additional anisotropics 
etc., should be made the subject of further work. In 
particular, the question of competing surface-couplings 
and partially compensated surfaces in finite antiferro- 
magnetic stacks should be addressed. In such systems, a 
competition between a genuine inhomogeneous spin-flop 
phase and twisted states takes place. 

In systems with intermediate anisotropics comparable 
to the indirect interlayer exchange within antiferromag- 
netic superlattices, one has to expect very complicated 
phase-diagrams. Still, such situations can be analysed by 
the micromagnetic methods developed here. However, it 
is vital to use clear concepts of magnetic phase transitions 
and clean definitions that designate the driving forces be- 
hind the varieties of field-driven reorientation processes 
in conflned antiferromagnets. We emphasize that the no- 
tion of a "surface spin-flop" is erroneous because the rel- 
evant magnetic energy terms that drive both the cant- 
ing instabilities at surfaces and the transition into the 
inhomogeneous spin-flop phase are not related to a bal- 
ance between effective anisotropies and Zeeman energy 
in these finite antiferromagnets. The transitions experi- 
enced by the type of finite antiferromagnets with non- 
compensated surfaces, as investigated here, are always 
related to the exchange cut. Finally, for the artificial 
antiferromagnetic systems composed of mesoscale ferro- 
magnetic units, the first-order transitions are of crucial 
importance. The phase-coexistence between states with 
finite magnetization will give rise to stable domain struc- 
tures and hysterctic behaviour in these systems owing to 
demagnetization effects. 
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